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ABSTRACT 

A simple model for the dynamics of dunes associated with vegetation is pro- 
posed. Using the model, formation processes of transverse dunes, parabolic dunes 
and elongated parabolic dunes according to two environmental factors: i)the 
amount of sand at the source, ii)the wind force, are simulated. The results have 
qualitative correspondence to the real counterpart, and the simplicity of the al- 
gorithm and the consequent easiness of the handling of this model provide us 
with wide applicability for the investigation of the complex interplay between 
vegetation and dunes. 

INTRODUCTION 

It is well known that the shape of dunes varies depending on several envi- 
ronmental factors surrounding them, like, the amount of available sand in each 
desert area, the wind directional variability, the vegetational condition covering 
sand surface, etc. Among all, the interplay between the dynamics of dunes and 
the growth of plants remains many aspects unsolved because of the complex na- 
ture of the dynamics of the system and the difficulty of the inspection of the 
theoretical hyposesis through direct observations. Still, many valuable observa- 
tional studies have been accumulated which provided us detailed environmental 
conditions at individual arid areas in which typical types of vegetated dunes are 
seen|T|, |^, |^, ^, |^, ^. One of the pioneering studies among them was made by 
Hack[0] in which he introduced a phase diagram to show the relation between 
the wind condition, available sand, the ratio of surface covered by plants, and 
the dominant type of dunes observed in each arid (or semiarid) area. Here, 
taking previous observational studies into consideration, we propose a minimal 
model which realizes the quahtative dynamics of dunes in mildly vegetated arid 
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areas. Thereafter, using the model, the formation processes of transverse dunes, 
parabohc dunes and elongate parabolic dunes according to environmental factors 
are simulated, that is, the Hack's phase diagram of vegetated dunes is numeri- 
cally testified. In the below we firstly make a brief explanation of our model and, 
after that, the outputs obtained through the simulation are shown. Finally we 
make a discussion on the meaning of the model. 

MODEL 



The model is a 2-dimensional(2D) lattice model, wherein two field variables,i)the 
local height h{i,j,n) of sand bed, and ii)the local density of vegetation c{i,j,n), 
are allocated at each site of the horizontally extending x N{N = 100) lattice. 
Here < i,j < N} and {n\0 < n} indicate, respectively, horizontal posi- 

tion and time step. Assume wind is constantly blowing in positive i direction, 
and each site covers a sufficiently wider area than that occupied by indi- 
vidual grains. Likewise is the unit time step which reflects a much longer time 
than individual saltation processes of sand. The above field variables are set to 
interact each other through the suppression factor apii^j^n) as explained below. 
Although many other quantities are considered to contribute to the whole dynam- 
ics of the system, we concentrate ourselves on extracting a simple set of relevant 
factors to investigate the essential dynamics of the whole system. Such method- 
ology for modeling the complex dynamics of dunes with a set of simple rules have 
recently been introduced by Warner and also by one of present authors[||, with 
which they succeeded to numerically reproduce the typical types of unvegetated 
dunes under conditions qualitatively corresponding to the observations made by 
Wasson and Hyde[0]. More previously, a pioneering model of the formation of 



vegetated dunes using a minimal set of rules has been introduced by Castro HlO 



Although his model is not applicable for the dynamics of dunes without vegeta- 
tion, neither, does reproduce a wide variety of morphology of dunes in vegetated 
areas, it successfully simulated a systematic change of the inter-dunes length of 
vegetated transverse dunes, also the change of their profiles. 

Present model is a hybrid model to investigate a wide class of morphodynamics 
of dunes with vegetation also without vegetation. Here, dynamical rules of sand 
movement are based on a previous work of the author 0, whereas the rules for 
time evolution of plants have many similarities to the model by Castro |T0[. 
I. For the evolution of c(z, j, ra), the local density of plants, we adopt a set of simple 
rules. In an extreme case where the height of sand surface remains unchanged 
with time, wherein plants can grow thick up to the saturation density without 
being cut away by the drastic defiation of ground or being buried by the rapid 
accumulation of sands, c(z, j, n) is allowed to increase linearly until the maximum 
value Cmax- On the other hand, if the temporal change of surface height is too 
fast, the growth rate of plants is suppressed or some of them may wither up and 
die, then c{i,j,n) decrease down to the minimum limit Cmin- To refiect the above 
situation, we use a discrete dynamics which is a sectional linear map as shown in 
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fig.l. Specifically, the dynamics is expressed as, 

C{i,j,n+1) = A{c{i,j,n) -b{i,j,n)) + Cmim {b < C < {Cmax - Cmin) / ^ + b) 

c{i,j,n+l) = Cmim {c<b) 

here b = b{i,j,n) = \h{i,j,n) — h{i,j,n — 1)|, c = c{i,j,n) and A is a constant 
to determine the growth rate of plants. 

II. for the evolution of h{i,j,n), discretized conservation law of h{i,j,n), 

h{i,j,n + 1) - h{ij,n) = Qinii,j,n) - Qoutihj,n) (1) 

holds, where Qinihjjn) is the total mass of sand coming into site at time 
step n, while Qoutih is the same quality leaving from at n. Both of 
the saltation flux and the creep flux contribute to Qin{i,j,n) and Qoutih 
Specifically the saltation flux caused by the grains leaving from at n is 

expressed by the production of its mass m and length / like, 

(lsai{i,j,n) = qo(tanh(yh{i,j,n)) + 1 .0) (tanh{-V h{i , j , n)) + 1.0 + a). (2) 

where Vh{i,j, n) means h{i,j, n) — h{i — n) and a is a constant to determine 
the bed-load in the windward slope of dunes. The above reflects the qualitative 
nature of wind, also the resulting saltation flux around dunes which is intensive 
in the windward particularly around the crest, whereas almost no flux in the 
lee side[|, |ll[(flg.l(b)). On the other hand, the flux by creep qcreepih is 
set proportional to the local gradient of sand surface. Although this may be a 
very crude approximation, it has some sense for a qualitative description of the 
morphodynamics of dunes as discussed afterward. The crucial effect caused by 
permitting the growth of plants is such that the sand transport sharply decreases 
when the cover ratio of sand surface by plants exceeds a critical value. To realize 
the situation in a simple expression, the suppression factor a^(i, j, n) is introduced 
like, 

ap{iJ,n) = l+(3{tanh{c-Cer)-l). (0 < /3 < 1/2) (3) 

Here, Ccr is the critical vegetation density over which the movement of sand 
sharply decreases, and f3 determines the maximum efficiency of suppression(flg.l(c)) 
The value of (3 depends on whether it is for the saltation flux or for the creep 
flux. With the suppression factor, the saltation/creep flux is forced to decrease 
as q{i,j,n) = a[^{i, j,n)'^q'{i, j,n) where q{i,j,n) and q'{i,j,n) are, respectively, 
the local flux by salt at ion/ creep with vegetation and without vegetation,. 

RESULTS 

Intending to compare the simulation outputs with the diagram by Hack, two 
kind of quantities are chosen as the control parameters; One is i)the amount of 
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sand at the source. Namely, the average height, < hgource >, of the sand surface 
at the source sites, = 1,2 1 < j < A^}. Specifically, uniformly random 

numbers between 2< hsource > and are allocated on these sites at each time 
step. The other control parameter is ii)the wind strength which should be a 
monotonically increasing function of saltation flux. Specifically, the variable go 
in eq.(2) is adopted as the index of wind strength. Note that the vegetation 
density, which is one of the axes in the diagram by Hack, is not adopted as a 
control parameter because it is rather a resultantly attained quantity after the 
above two control parameters are fixed. Simulations are initiated from fiat sand 
surface except the source area, namely, h{i,j, 0) = holds except the source sites, 
while the initial vegetation density c{i,j, 0) at each site is set random around the 
average value < c(z,j, 0) > which is between Cmax and Cmin- Note the boundary 
condition in j direction, which is perpendicular to the wind direction, is set 
periodic and the leeward boundary in i direction is set as the free boundary. 
Below the initial level of sand surface, erosion of sand surface is inhibited to 
realize the existence of the hard ground or the ground water table. Using these 
rules, spontaneous formation processes of typical types of dunes are observed. 

The results are; First of all, when the wind force is too weak, regardless of 
the amount of sand supply at the source, clearly shaped dune will not appear in 
the system(fig.2(a)). It is also the cases where the amount of sand supply is too 
small. On the other hand, with more than a certain strength of wind force and 
a certain amount of sand supply, two types of clearly shaped active dunes are 
formed; 

i) When sufficiently large amount of sand is supplied under rather strong wind, 
transverse dunes, barchan dunes or both of them will dominate in the system. In 
more detail, small parabolic dunes formed just lee of the sand source soon develop 
into barchan dunes, which connect each other to grow larger as they move, to 
form transverse dunes the crests of them extend, roughly, perpendicular to the 
direction of wind(fig.2(b)). 

ii) If the amount of supplied sand or the wind force is slightly less than the above, 
parabolic dunes will prevail in the system(fig.2(c)(d)). They have arms extending 
to the windward direction. At the center of each parabolic dune, hollow is devel- 
oped in which surface erosion proceeds down to the level of hard ground to serves 
the dune's nose with sand. There we can see the inclination that the length of 
the arms varies depending on the wind force, namely, the arms extend the longer 
under the stronger wind. Especially in the case with rather small amount of sand, 
elongate parabolic dunes[§[ with the shapes like hair-pins will appear(fig.2(d)). 
The arms of such elongate dunes, if without their noses, seem like pairs of linear 
dunes. 

DISCUSSION 

The above results are summarized in a phase diagram as shown in fig. 4. This 
diagram is not directly compared with that of Hack because of the less number 
of axes in our diagram. Also, in the diagram, ambiguous area is left in which 
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condition irregular mounds of sand are formed, which are not easily categorized 
into other dominant types of dunes. Still, it has qualitatively good correspondence 
to the previous observational studies including that of Hack. Especially, the 
systematic change of the dunes' morphology from the transverse(or barchan) 
dunes to the parabolic or the elongate parabolic dunes according to the amount 
of available sand (or to the wind force), is clearly simulated. Moreover no dune 
formation is realized under too weak wind nor is under too small amount of sand 
supply. The results indicate that this simple model contains intrinsic dynamics 
by which the morphodynamics of the vegetated dunes is decisively affected, and 
that through this model theoretical hypotheses on the formation dynamics of 
vegetated dunes can be testified. 

On the other hand, present study leaves several, possibly crucial, situations 
untreated. Firstly, the wind is always blowing in one direction. This condition 
seems very hard for the spontaneous growth of linear dunes considering the cases 
of unvegetated dunes. Consequently, straightly extending sand dunes/ridges ap- 
pear only in behind parabolic dunes as the traces of their arms. This is one of 
the possible origin of vegetated linear dunes But it seems not the general sce- 
nario for the formation of linear dunes, thus, the study under complicated wind 
regimes is one the next issues. Secondly, the existence of the angle of repose was 
not explicitly incorporated in the model. Here, instead, all the dynamics along 
sand surface was modeled by the diffusion-like creep process. Of course such 
a rough simplification can cause the imperfect reproduction of the real counter- 
parts, like inhomogeneous angle of slope at the lee face of a dune unlike the actual 
slip face. On the problem, in the previous studied, we demonstrated that even 
such a simplified modeling can work effectively, at least for the investigation of 
macroscopic morphodynamics of unvegetated dunes, like, simulating distinguish- 
able barchan, self, star and other types of dunes. Also in the present simulations 
with vegetation, distinguishable parabolic and other types of dunes successfully 
appeared according to the corresponding situations in the semiarid desert areas. 
However, for the investigation of more quantitative aspects of the system, the 
explicit introduction of the angle of repose into the model is required, and it also 
be an important step to the inclusive understanding of the system. 

Notwithstanding all, we believe the present study using such a minimal model 
is one of the effective way to understand the complex interplay between dunes 
and vegetation. Systematic search of complex systems using such bold simplifi- 
cations would , more or less work as the complementary way which fill the big 
blank between the observational studies of such systems and the, otherwise hardly 
testified, theoretical hyposesis. 
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Figure Captions 



Fig.l Schematic explanations of the present model; (a-I)The Map c{i,j,n) — 
c{i,j,n + 1) to describe the discretized time evolution of plants density. Without 
the temporal change b{i,j,n) of surface height, plants density monotonically in- 
creases up to the saturation value Cmax- (a-II)With more than a certain speed of 
surface rise or deflation, plants at the surface are, more or less, damaged because 
they are buried or cut away, then, c{i,j, n) will decrease with time with the lower 
limit Cmin- b)Local sand flux by saltation qsai of eq.(2) is given as a function of 
local slope Vh of sand surface in the direction of wind, which reflects the ob- 
servational fact, namely, large bed load in the windward particularly around the 



crest (indicated by the arrow), whereas sharply it decreases in the lee side[|], |TT 
c) Suppression factor a of sand flux as a function of local density of vegetation c. 
Above a critical density c^r of vegetation, the flux of sand drastically decreases 
as described in eq.(3). 

Fig. 2 Snapshots of simulated dunes under various pairs of control parameters, 
wind force and sand supply. The left part in each figure shows the spatial distri- 
bution of vegetation density d{i, j,n). The darker tone indicates the more densely 
vegetated place, whereas white parts indicates the areas with bare sand surfaces. 
The right part in each figure shows the surface height distribution h{i,j.n), where 
the darker position means the higher surface. Sand is supplied from the most 
windward 2 rows. Steady wind is blowing from the left to the right, (a) When the 
wind force is too weak no distinguishable dune is formed. It is also the case for 
too small amount of sand supply at the source. (b)Under rather strong wind with 
sufficient amount of sand supply transverse dunes prevail while small parabolic 
dunes are seen just lee of the sand source. The latter will soon grow up to the 
former, (c) Parabolic dunes, the arms of which extend in the windward direction, 
are formed under mildly blowing wind with intermediate amount of sand at the 
source. (d)If the amount of sand is comparatively small within this regime, thin 
and long parabolic dunes, namely, elongate parabolic dunes will grow. They look 
like rather linear dunes if without their noses. 



Fig. 3 Phase diagram to show the dominant types of dunes under various pairs 
of control parameters. The alphabets in the diagram indicate the conditions 
corresponding to respective snapshots in fig.2. Symbols As and "s indicate, 
respectively, parabolic dunes, and, barchan(or transverse) dunes, whereas xs 
mean the conditions for no dune formation. The As accompanied by * mean the 
conditions for the development of rather thin parabolic dunes with long arms. At 
the condition with the symble □, many irregular mounds are formed which are 
not clearly categorized as particular type of dunes. 
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